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Abstract 

Presented  in  this  paper  is  a  computational  analysis  of  a  membraneless  microfluidic  fuel  cell  that  uses  the  laminar  nature  of  microflows  to 
maintain  the  separation  of  fuel  and  oxidant  streams.  The  fuel  cell  consists  of  a  T-shaped  microfluidic  channel  with  liquid  fuel  and  oxidant 
entering  at  separate  inlets  and  flowing  in  parallel  without  turbulent  or  convective  mixing.  Recent  experimental  studies  have  established 
proof-of-concept  of  such  fuel  cells  and  have  also  shown  that  their  performance  is  greatly  limited  by  poor  fuel  utilization.  Improving  fuel 
utilization  while  minimizing  fuel-oxidant  mixing  in  microfluidic  fuel  cells  is  the  focus  of  this  study.  A  concise  electrochemical  model  of 
the  key  reactions  and  appropriate  boundary  conditions  are  presented  in  conjunction  with  the  development  of  a  computational  fluid  dynamic 
(CFD)  model  of  this  system  that  accounts  for  coupled  flow,  species  transport  and  reaction  kinetics.  3D  numerical  simulations  show  that  the 
fuel  cell  is  diffusion  limited,  and  both  microchannel  and  electrode  geometry  play  key  roles  in  system  performance.  Three  cross-sectional 
geometries  are  investigated,  and  a  high  aspect  ratio  rectangular  geometry  results  in  a  two-fold  increase  in  fuel  utilization  compared  to  a  square 
geometry  with  the  same  hydraulic  diameter.  By  tailoring  the  flow  rate  to  the  axial  length  of  the  fuel  cell,  fuel  utilization  is  increased  to  23%. 
Using  the  numerical  simulation  to  guide  the  electrode  design  process,  an  extended  tapered-electrode  design  is  proposed.  Simulations  of  the 
tapered-electrode  microfluidic  fuel  cell  demonstrate  a  fuel  utilization  of  over  50%. 

©  2004  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Micro-sized  power  sources  are  needed  for  small,  portable 
devices  capable  of  operating  for  long  periods  without 
recharging,  such  as  cell  phones  and  laptop  computers  through 
to  more  specialized  devices  such  as  remote  sensors,  global  po¬ 
sitioning  devices,  and  in  vivo  diagnostic  medical  devices.  It  is 
predicted  [1]  that  battery  technology  will  not  keep  pace  with 
these  growing  portable-power  demands,  particularly  with  the 
next  wave  of  wireless  technology,  broadband  mobile  comput¬ 
ing.  Microstructured  fuel  cells  have  the  potential  to  bridge  the 
gap  between  battery  technology  and  growing  portable-power 
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demands,  by  providing  longer  duration  per  weight  and  vol¬ 
ume,  quasi-instant  recharge,  and  constant  discharge.  Recent 
developments  in  the  field  have  benefited  from  micromanufac¬ 
turing  technology  and  biological  and  chemical  lab-on-a-chip 
concepts  that  have  been  a  driving  force  behind  the  recent  de¬ 
velopments  of  many  microfluidic  devices  [2] .  Other  emerging 
application  areas  for  microfluidics  include  electronic  cooling, 
aerospace,  and  printing  [3]. 

Since  all  reactions  taking  place  in  fuel  cells  are  surface 
based,  the  increase  in  the  surface-to-volume  ratio  accompa¬ 
nying  miniaturization  leads  to  a  fundamental  improvement  in 
power  density  [4,5].  There  are,  however,  mechanical  limits  to 
the  miniaturization  of  conventional  fuel  cells  [6] :  machining 
graphite  bipolar  plates  becomes  difficult,  and  decreasing  the 
size  of  the  membrane  and  substrate  decreases  their  support- 
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Nomenclature 

A  pre-exponential  factor 

C\  concentration  of  cell  i 

C*mxed  concentration  of  perfectly  mixed  solution  at 
cell  i 

Cf  concentration  of  solution  at  cell  i  in  the  absence 
of  mixing 

Di  diffusion  coefficient 

E  potential  at  electrode 

Ea  activation  energy 

Efev  reversible  potential 

F  Faraday  constant 

h  channel  height 

Ji  flux  of  species  i 

k  rate  constant 

n  number  of  electrons 

N  total  number  of  cells 

P  static  pressure 

R  gas  constant 

RCq\\  ionic  resistance  of  electrolyte 

Ri  net  rate  of  production  of  species  i 

t  time 

T  temperature 

U  average  fluid  velocity 

v  rate  law 

v  velocity  vector 

Av  diffusive  mixing  region  width 

y  distance  fluid  travels  downstream 

Yi  local  mass  fraction  of  species  i 

Greek 

a  charge  transfer  coefficient 

r]  overpotential 

p  fluid  density 

v  kinematic  viscosity 

X  apparent  reaction  order 

Chemical  formulas 
CO2  carbon  dioxide 

(COOH)ads  adsorbed  COOH 
e_  electron 

H+  proton 

HCOOH  formic  acid 
O2  oxygen 

Pt  platinum 

Pt°  occupied  platinum  site 


ing  strength.  In  addition,  the  performance  of  miniaturized 
conventional  fuel  cells  is  limited  by  the  Ohmic  polarization 
introduced  by  the  membrane  and  related  water  management 
issues.  Novel  architectures  are  crucial  for  the  success  of  mi¬ 
croscale  fuel  cells  [4] . 


Recently  introduced  membraneless  microfluidic  fuel  cells 
[7-9]  take  advantage  of  the  laminar  nature  of  microflows  to 
maintain  the  separation  of  fuel  and  oxidizer  streams  without 
the  use  of  a  membrane.  A  schematic  of  a  microfluidic  fuel 
cell  and  its  operation  are  given  in  Fig.  1 .  The  geometry  is  that 
of  a  T-intersection,  or  T-mixer,  which  is  commonly  used  in 
analytical  microfluidic  chips.  The  cross-stream  mixing  rate 
in  such  T-intersections  is  diffusion  limited,  and  many  studies 
have  focused  on  increasing  the  mixing  rate  in  such  geometries 
[2] .  The  microfluidic  fuel  cell,  on  the  other  hand,  exploits  the 
nature  of  this  flow  to  achieve  the  separation  of  fuel  and  oxi¬ 
dant  streams.  The  separation  of  fuel  and  oxidant  is  required 
to  restrict  reactions  of  oxidation  and  reduction  to  the  appro¬ 
priate  electrode.  Fuel  is  introduced  at  one  inlet,  and  oxidant  is 
introduced  at  the  second  inlet.  Electrodes  are  placed  along  the 
walls  to  complete  the  fuel  cell,  and  oxidation  at  the  anode  and 
reduction  at  the  cathode,  together  provide  the  cell  potential. 
The  fluid  facilitates  protonic  conduction  from  one  electrode 
to  the  other,  and  the  electrons  generated  at  the  anode  take 
an  external  path  through  an  applied  load.  The  length  of  the 
reaction  channel  is  limited  by  the  mixing  of  the  two  streams. 
The  operation  of  the  fuel  cell  will  begin  to  fail  when  the  two 
streams  become  mixed  to  the  point  that  oxidation  and  reduc¬ 
tion  are  no  longer  restricted  to  the  appropriate  electrodes. 

Choban  et  al.  [9]  were  the  first  to  demonstrate  a  mem¬ 
braneless  fuel  cell  using  formic  acid  and  oxygen  as  reac¬ 
tants.  They  demonstrated  that  when  two  streams  are  flowing 
in  parallel  in  the  laminar  regime,  the  streams  remain  sepa¬ 
rated,  eliminating  the  need  for  a  membrane.  Ferrigno  et  al. 
[8]  demonstrated  a  millimeter-scale  fuel  cell  using  vanadium 
as  reactants.  The  advantage  of  this  design  is  that  it  used  the 
same  species  (vanadium)  as  fuel  and  oxidant,  which  can  be 
regenerated  from  a  mixture  of  the  products  [8].  Choban  et 
al.  [7]  reported  a  Y-shaped  microfluidic  membraneless  fuel 
cell  system  using  formic  acid  and  oxygen  as  reactants,  which 
reached  a  current  density  of  0.4  mA  cm-2.  The  advantage  of 
this  system  is  that  it  uses  acidic  solutions  to  minimize  the 
protonic  resistance  in  the  fluid.  Both  systems  [7,8]  were  re¬ 
ported  to  be  diffusion  limited.  In  this  work,  a  microchannel 
geometry  with  a  high  aspect  ratio  in  the  cross- stream  direc¬ 
tion,  similar  to  that  of  Ferrigno  et  al.  [8],  is  investigated, 
as  well  as  a  square  geometry  similar  to  that  of  Choban  et 
al.  [7].  Microfluidic  fuel  cells  have  several  advantages  over 
conventional  fuel  cells:  eliminating  the  membrane  removes 
related  Ohmic  losses,  water  management  and  sealing  issues; 
and  since  the  fuel  and  oxidant  streams  flow  together  in  the 
same  channel  network,  the  fuel  cell  size  is  reduced,  the  de¬ 
sign  is  simplified,  and  the  manifolding  requirements  are  also 
reduced.  In  addition,  research  in  this  area  can  capitalize  on 
recent  advances  in  numerical,  experimental,  and  microfabri¬ 
cation  techniques  [10].  Microfluidic  fuel  cells  have  the  fur¬ 
ther  advantage  of  using  liquid  fuels,  which  have  high  energy 
densities  compared  to  gaseous  fuels  [1],  particularly  impor¬ 
tant  in  the  context  of  portable-power  applications.  However, 
current  designs  show  relatively  poor  fuel  utilization,  on  the 
order  of  1%  [7]  to  10%  [8].  MicroChannel  geometries  and 
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Fig.  1.  Numerical  simulation  and  schematic  of  a  microfluidic  fuel  cell  with  fuel  and  oxidant  entering  at  separate  inlets  and  remaining  separated  as  they  flow 
downstream.  The  flow  rates  and  length  scales  were  selected  to  illustrate  the  mixing  process. 


tailored  electrode  placements  are  investigated  to  consider¬ 
ably  improve  current  designs.  Modelling  can  provide  insight 
into  the  parameters  and  geometry  required  to  improve  the 
performance  of  this  technology. 

In  this  work  a  numerical  analysis  of  a  membraneless  mi¬ 
crofluidic  fuel  cell  is  presented.  This  work  is  one  of  the  first 
computational  fluid  dynamic  (CFD)  analyses  of  this  technol¬ 
ogy.  The  model  accounts  fully  for  three-dimensional  con¬ 
vective  transport  in  conjunction  with  anodic  and  cathodic 
reaction  kinetics.  Appropriate  boundary  conditions  for  the 
CFD  modelling  of  this  system  are  developed  here  and  im¬ 
plemented  into  the  numerical  model.  A  numerical  investiga¬ 
tion  of  the  coupled  flow,  species  transport  and  electrochem¬ 
ical  aspects  in  this  system  is  conducted.  The  results  provide 
insight  into  the  running  parameters  and  both  microchannel 
and  electrode  geometries  required  to  achieve  significantly 
improved  performance.  Finally,  using  the  numerical  simu¬ 
lation  to  guide  the  electrode  design  process,  an  extended 
tapered-electrode  design  is  proposed,  and  its  performance  is 
investigated. 


2.  Hydrodynamic  and  mass  transport  model 

A  three-dimensional  CFD  model  is  applied  to  simulate  the 
coupled  flow,  species  transport,  and  chemical  aspects  of  the 
microfluidic  fuel  cell.  Due  to  the  moderate  liquid  velocities, 
the  internal  heating  due  to  viscous  dissipation  is  ignored,  and 
an  isothermal  system  is  assumed  [3] .  The  following  equations 
are  solved  for  laminar  flow  in  an  inertial  reference  frame 
at  steady  state  for  incompressible  and  isothermal  fluid  flow. 
Neglecting  body  forces,  the  Navier  Stokes  and  continuity 
equations  for  incompressible  and  isothermal  flow  are  given 
by  [11]: 

dl)  1 

- F  V  •  (vv)  —  — VP  +  V  •  (uVu)  (1) 

dt  p 

V  •  v  =  0  (2) 


where  P  is  the  static  pressure,  p  the  fluid  density,  v  the  velocity 
vector,  and  v  the  kinematic  viscosity.  The  conservation  of 
species  equation  is  given  by 

V  •  (pvYi)  =  -V  •  %  +  Ri  (3) 

where  Fz  is  the  local  mass  fraction  of  species  i  and  Ri  is  the  net 
rate  of  production  of  species  i  by  chemical  reaction.  For  the 
dilute  approximation  that  is  used  in  this  model,  the  diffusion 
flux  of  species  i  is  given  by  Fick’s  law: 

Ji  =  -pDiVYi  (4) 

where  Df  is  the  diffusion  coefficient  for  species  i  into  the 
mixture. 

3.  Reaction  model 

3.1.  Formic  acid  oxidation 

Formic  acid  is  an  attractive  fuel  for  fuel  cells,  as  it  has  been 
reported  to  be  electrochemically  more  active  than  methanol 
[12].  The  oxidation  of  formic  acid  is  a  well-studied,  yet  com¬ 
plex  chemical  reaction.  The  challenge  is  to  provide  a  reason¬ 
able  model  for  this  highly  non-linear,  dual-pathway  mech¬ 
anism  [13,14].  A  platinum  (Pt)  catalyst  is  modelled  due  to 
its  well- studied  nature  in  regards  to  formic  acid  oxidation  as 
well  as  use  in  conventional  fuel  cells. 

The  oxidation  of  formic  acid  (HCOOH)  on  Pt  occurs  via  a 
dual  pathway  mechanism  [15].  With  the  proper  catalyst,  the 
oxidation  of  formic  acid  favors  the  direct  pathway,  the  dehy¬ 
drogenation  reaction  [15-18]  through  which  carbon  dioxide 
(CO2)  is  formed  directly  according  to  the  generally  accepted 
mechanism  proposed  by  Capon  and  Parsons  [19]: 

Pt°  +  HCOOH  -*  (COOH)ads  +  H+  +  e“  (5) 

(COOH)ads  -*  C02  +  H+  +  e"  +  Pt°  (6) 

The  rate  determining  step  occurs  between  HCOOH  and  the 
Pt  electrode  with  the  transfer  of  one  electron,  as  shown  in  Eq. 
(5)  [14,16].  Modelling  the  oxidation  of  formic  acid  entails  the 
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formulation  of  the  rate  constant  and  rate  law  governing  the 
reaction.  The  rate  constant  in  a  chemical  reaction  is  given  by 
the  Arrhenius  equation: 

k  =  Ae~E*IRT  (7) 

where  A  is  the  pre-exponential  factor,  Fa  the  activation  en¬ 
ergy,  R  the  gas  constant,  and  T  the  temperature.  The  chemical 
reaction  is  also  governed  by  the  rate  law,  which  is  unique  to 
each  chemical  reaction  [20] .  The  rate  law  can  also  be  affected 
by  the  coverage  of  platinum  sites  by  carbon  monoxide  and 
by  (COOH)ads  from  the  second  step  of  the  direct  pathway, 
Eq.  (6).  However,  since  the  second  step  in  the  dehydrogena¬ 
tion  reaction  is  believed  to  be  faster  than  the  rate  determining 
step  [14],  it  is  assumed  that  the  coverage  of  platinum  sites 
by  (COOH)ads  is  negligible.  It  is  also  assumed  that  the  fuel 
cell  will  be  operated  at  a  potential  where  poisoning  is  not 
an  issue,  so  that  coverage  by  carbon  monoxide  is  negligible. 
Therefore,  the  rate  law  is  given  by  [14]: 

v  =  k  [HCOOH] x  exp  f 

where  [HCOOH]  is  the  concentration  of  HCOOH  in  the  so¬ 
lution,  x  the  apparent  reaction  order,  a  the  charge  transfer 
coefficient,  n  the  number  of  electrons  released  at  that  partic¬ 
ular  step,  F  the  Faraday  constant,  and  rj  the  overpotential. 

In  the  rate  determining  step,  one  electron  is  released,  there¬ 
fore  ft  —  1  in  Eq.  (8).  The  reversible  potential  of  formic  acid 
oxidation  is  EYQV  =  —0.07  V  (standard  hydrogen  electrode), 
so  for  £  =  0.6  V  reported  by  Fonseca  et  al.  [16],  the  overpo¬ 
tential  is  given  by  r\-  0.67  V.  According  to  the  law  of  mass 
action  [21],  the  rate  determining  step  is  a  first-order  reac¬ 
tion,  making  /=  1.  According  to  Capon  and  Parsons  [19], 
the  charge  transfer  coefficient  for  the  oxidation  of  formic 
acid  on  platinum  is  a  =  0.51.  Similarly,  Jiang  and  Kucernak 
[14]  reported  the  experimentally  determined  Tafel  slope  to 
be  132  mV  dec-1  at  60  °C,  which  resulted  in  a  charge  trans¬ 
fer  coefficient  of  a  =  0.497.  The  charge  transfer  coefficient 
assumed  for  our  model,  a  =  0.497,  is  in  keeping  with  these 
works  [14,19]. 

For  the  rate  determining  step  of  the  formic  acid  oxi¬ 
dation  reaction,  Fonseca  et  al.  [16]  reported  a  rate  con¬ 
stant  of  k=10~2cms_1  on  a  platinum  catalyst  modified 
with  lead  adatoms  at  an  electrode  potential  of  0.6  V  (nor¬ 
mal  hydrogen  electrode)  at  18±2°C.  Pletcher  and  Solis 
[17]  reported  a  rate  constant  for  the  rate  determining  step 
ofk=  1.0  ±0.5  x  10  2  cm  s  1  at  a  potential  of  0.6  V  at  room 
temperature  for  the  oxidation  of  formic  acid  on  a  platinum  an¬ 
ode  with  lead  adatoms.  Pletcher  and  Solis  [17]  and  Fonseca  et 
al.  [16]  both  used  the  same  coverage  of  50%  of  lead  adatoms 
on  their  Pt  surfaces,  and  obtained  rate  constants  within  the 
same  order  of  magnitude.  The  rate  constant  assumed  for  our 
model,  k-  lx  10-2  cms-1,  is  in  keeping  with  these  works 
[16,17]. 


3.2.  Oxygen  reduction 

The  protons  released  from  the  oxidation  reaction  at  the  an¬ 
ode  diffuse  through  the  solution  in  the  cross-stream  direction 
towards  the  cathode.  At  the  cathode,  the  oxygen  (O2)  from 
the  oxidant  stream,  the  protons  and  the  electrons  that  travel 
through  the  external  circuit  combine  to  produce  water  in  the 
oxygen  reduction  reaction: 

02  +  4H+  +  4e”  2H20  (9) 

Neglecting  the  structural  sensitivity  of  the  Pt  surface,  a  rate 
constant  of  k  =  250cms-1  is  assumed  for  our  model  [22]. 
Based  on  this  assumption,  poisoning  of  the  platinum  surface 
is  neglected,  and  the  rate  law  can  be  written  as 

V  =  k[02\  exp  (10) 

where  [O2]  is  the  concentration  of  O2  in  the  solution.  The 
reversible  potential  for  oxygen  reduction  at  standard  state  is 
1.46  V.  From  the  definition  of  overpotential,  r\-  —0.66  V  at 
an  electrode  potential  of  0.8  V.  At  this  electrode  potential,  the 
charge  transfer  coefficient  is  a  =  0.5  [22]. 

Both  the  formic  acid  oxidation  and  oxygen  reduction  re¬ 
actions  were  assumed  to  be  first-order  reactions.  Several  as¬ 
sumptions  have  been  made  in  order  to  present  a  reasonable 
model  for  the  implementation  of  redox  reactions  into  the  com¬ 
putational  model.  Surface  poisoning  has  been  neglected,  the 
dominance  of  rate  determining  steps  has  been  assumed,  and 
electrochemical  parameters  have  been  approximated  based 
on  empirical  results. 

4.  Membraneless  microfluidic  fuel  cell 

The  model  was  implemented  in  Fluent,  a  commercial 
finite-volume  based  CFD  package.  Dual  processors  were 
used  to  calculate  the  solutions  generated  from  a  hexagonal 
mesh,  with  typical  run  times  of  2h.  In  all  simulations,  a  hy¬ 
draulic  diameter  of  100  p,m  was  used  for  the  inlet  and  outlet 
channels.  Simulations  involving  the  hydrodynamic  and  mass 
transport  model  are  presented,  followed  by  simulations  in¬ 
cluding  the  reaction  model. 

The  fuel  and  oxidant  cross-stream  diffusive  mixing  was 
quantified  with  the  hydrodynamic  and  mass  transport  model. 
The  operation  of  the  fuel  cell  will  begin  to  fail  when  the  two 
streams  become  mixed  to  the  point  that  oxidation  and  re¬ 
duction  are  no  longer  restricted  to  the  appropriate  electrodes. 
The  extent  to  which  the  fuel  and  oxidant  become  diffusively 
mixed  is  proportional  to  the  diffusivity  of  the  fuel  in  the  sol¬ 
vent.  The  diffusivity  also  affects  the  mass  transport  of  the  fuel 
to  the  electrodes.  The  binary  mixture  is  assumed  dilute,  and 
the  diffusion  coefficient  is  constant  throughout.  A  diffusion 
coefficient  of  D  =  5  x  10  10m2s  1  is  assumed,  which  is  a 
typical  value  for  the  diffusion  of  relatively  small  molecules 
in  an  aqueous  solution  [2,20]. 


A.  Bazylak  et  al.  /  Journal  of  Power  Sources  143  (2005)  57-66 


61 


Fig.  2.  Liquid  transport  in  the  microfluidic  fuel  cell:  (a)  schematic  illustrating  the  location  of  the  X-Z  cross  section  of  interest  for  the  axial  velocity  contour 
plot  of  the  square  geometry  (100  |xm  x  100  pun),  (b)  corresponding  axial  velocity  contour  plot,  (c)  schematic  illustrating  the  location  of  the  X-Z  cross  section 
of  interest  for  the  axial  velocity  contour  plot  of  rectangular  geometry  (550  |xm  x  55  pun),  and  (d)  corresponding  axial  velocity  contour  plot. 


Pressure  was  applied  to  drive  the  flow,  and  a  no-slip  ve¬ 
locity  condition  is  applied  to  the  walls.  These  boundary  con¬ 
ditions  result  in  the  formation  of  boundary  layers  with  steep 
cross-stream  velocity  gradients,  apparent  in  the  velocity  pro¬ 
files  plotted  in  Fig.  2.  These  cross-stream  velocity  gradients 
greatly  affect  the  diffusive  mixing.  Near  the  walls,  the  fluid 
flows  relatively  slowly.  As  a  result,  cross-stream  mixing  is 
more  pronounced  in  the  near- wall  regions,  and  the  mixing 
region  takes  on  an  hour-glass  shape  as  shown  in  Fig.  1.  The 
calculated  width  of  the  diffusive  mixing  region  at  the  walls  of 
the  fuel  cell  is  in  agreement  with  the  theoretically  predicted 
[23]  and  experimentally  validated  [8]  trend  represented  by 


the  following  expression: 


(11) 


where  h  is  the  channel  height,  y  the  distance  the  fluid  flows 
downstream,  and  U  the  average  flow  speed.  Down  the  cen¬ 
ter  of  the  channel,  the  diffusive  mixing  region  grows  more 
slowly,  following  a  one-half  power  scaling  [23].  Several 
works  have  focused  on  the  inclusion  of  electrokinetic  phe¬ 
nomena  in  microfluidic  systems  [24-26],  but  due  to  the  high 
ionic  concentrations  employed  here,  electrokinetic  effects 


Geometry  3 


Fig.  3.  A  schematic  diagram  of  the  microfluidic  fuel  cell  system:  (a)  the  location  of  the  X-Z  cross  section  of  interest,  (b)  Geometry  1  is  square  with  electrodes 
placed  on  the  top  and  bottom  channel  surfaces,  (c)  Geometry  2  is  rectangular  with  electrodes  placed  on  the  top  and  bottom  channel  surfaces,  and  (d)  Geometry 
3  is  rectangular  with  electrodes  placed  along  the  left  and  right  walls  of  the  channel. 
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Cross  stream  direction  (pin) 


Fig.  4.  Formic  acid  bulk  concentration  normal  to  the  anode  wall  for  various  downstream  positions  (at  the  anode  midpoint).  The  concentration  is  reduced  to 
zero  at  the  electrode  surface,  indicating  that  the  reaction  kinetics  are  fast  compared  to  the  diffusion  in  this  system. 


are  negligible.  The  challenges  in  modelling  microfluidic  fuel 
cells  lie  in  capturing  the  three-dimensional  mixing  dynamics 
and  the  reaction  and  electrode  kinetics. 

Fig.  3  shows  a  schematic  of  the  microfluidic  fuel  cell  with 
different  geometries,  all  of  which  have  a  hydraulic  diameter 
of  100  jam.  Geometry  1  shown  in  Fig.  3(b)  and  Geometry  2 
shown  in  Fig.  3(c)  have  electrodes  placed  along  the  top  and 
bottom  walls  with  a  20  (am  spacing  between  the  electrodes  to 
account  for  the  diffusive  mixing  regions.  A  high  aspect  ratio 
is  used  for  Geometry  2,  and  a  low  aspect  ratio  is  used  for 
Geometry  3,  whose  electrodes  are  placed  along  the  left  and 
right  walls  of  the  channel. 

Fig.  4  shows  the  concentration  profiles  of  formic  acid  at 
the  midpoint  of  the  anode  when  the  electrode  kinetics  were 
included  in  the  model.  The  concentration  was  reduced  to  zero 
at  the  electrode  surface,  indicating  that  the  reaction  kinetics 
were  fast  compared  to  the  diffusion  in  this  system.  This  result 
is  in  keeping  with  experimental  results  [7,8].  Since  the  fuel 
cell  is  diffusion  limited,  the  geometry  of  the  microchannel 
plays  a  predominant  role  in  the  efficiency  of  the  cell. 

Although  the  fuel  and  oxidant  streams  do  not  experience 
turbulent  or  convective  mixing,  they  do  experience  diffusive 


mixing.  The  formic  acid  and  oxygen  concentration  bound¬ 
ary  layers  developed  at  the  electrode  surfaces,  along  with  the 
diffusive  mixing  regions  are  illustrated  in  Figs.  5-7,  for  Ge¬ 
ometries  1-3,  respectively.  Figs.  5-7 (b)  shows  the  formic  acid 
mass  fraction  contour  plot,  showing  the  diffusive  mixing  be¬ 
tween  fuel  and  oxidant  streams  at  the  outlet  in  the  absence  of 
electro-oxidation  reactions.  As  expected,  the  diffusive  mix¬ 
ing  region  takes  on  an  hour-glass  shape  due  to  the  slower 
velocities  at  the  walls.  Figs.  5-7(c)  shows  the  formic  acid 
mass  fraction  contour  plots  combining  the  effect  of  diffusive 
mixing  and  the  electro-oxidation  at  the  anodes.  A  thicker 
concentration  boundary  layer  is  apparent  near  the  wall  ad¬ 
jacent  to  the  electrode  surface  and  is  due  to  the  combined 
wall/end  effects.  Similarly,  Figs.  5-7(d)  shows  the  oxygen 
mass  fraction  contour  plots  combining  diffusive  mixing  and 
the  electro-reduction  at  the  cathode.  Due  to  a  constant  dif¬ 
fusion  coefficient  and  diffusion  limited  electrodes,  the  de¬ 
pletion  and  mixing  of  oxygen  in  (d)  is  a  mirror  image  of 
the  depletion  and  mixing  of  formic  acid  in  (c).  In  each  case 
the  Reynolds  number  is  20,  based  on  the  hydraulic  diameter, 
which  is  100  ptm  in  each  case.  A  length  of  6  mm  downstream 
was  chosen  as  the  outlet  in  each  case  for  comparison  pur- 
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Fig.  5.  Depletion  and  mixing  of  HCOOH  and  O2  in  the  microfluidic  fuel  cell  with  Geometry  1 :  (a)  the  location  of  the  X-Z  cross  section  of  interest,  (b)  HCOOH 
mass  fraction  contour  plot  showing  diffusive  mixing  between  the  fuel  and  oxidant  streams  without  electro-oxidation  reactions,  (c)  HCOOH  mass  fraction 
contour  plot  showing  combined  diffusive  mixing  and  electro-oxidation  at  the  anodes  (top  and  bottom,  left),  and  (d)  O2  mass  fraction  contour  plot  showing 
combined  diffusive  mixing  and  reduction  at  the  cathodes  (top  and  bottom,  right). 
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Fig.  6.  Depletion  and  mixing  of  HCOOH  and  O2  in  the  microfluidic  fuel  cell  with  Geometry  2:  (a)  the  location  of  the  X-Z  cross  section  of  interest,  (b)  HCOOH 
mass  fraction  contour  plot  showing  diffusive  mixing  between  the  fuel  and  oxidant  streams  without  electro-oxidation  reactions,  (c)  HCOOH  mass  fraction 
contour  plot  showing  combined  diffusive  mixing  and  electro-oxidation  at  the  anodes  (top  and  bottom,  left),  and  (d)  O2  mass  fraction  contour  plot  showing 
combined  diffusive  mixing  and  reduction  at  the  cathodes  (top  and  bottom,  right). 


poses.  Fuel  utilization  was  determined  from  the  ratio  of  fuel 
mass  flow  rate  at  the  outlet  to  the  fuel  mass  flow  rate  at  the 
inlet.  With  an  inlet  velocity  of  0. 1  m  s-1 ,  the  square  geometry 
provides  a  fuel  utilization  of  3%,  whereas  the  rectangular  ge¬ 
ometries  with  aspect  ratios  of  1:10  provide  a  fuel  utilization 
of  8%.  These  results  are  in  agreement  with  recent  experimen¬ 
tally  determined  fuel  utilizations  on  the  order  of  1  %  for  the 
square  geometry  presented  by  Choban  et  al.  [7]  and  10%  for 
the  rectangular  geometry  presented  by  Ferrigno  et  al.  [8].  Al¬ 
though  Geometries  2  and  3  provide  similar  fuel  utilization, 
there  is  significantly  less  cross-stream  fuel-oxidant  mixing  in 
Geometry  2,  illustrated  in  Fig.  6(b).  Both  geometries  exhibit 
a  mixed  region  of  similar  width,  but  the  volume  of  fluid  in  the 
mixed  region  is  significantly  less  in  Geometry  2.  The  degree 
to  which  the  reactants  are  mixed  is  an  indication  of  how  much 
potential  there  is  for  improved  fuel  utilization.  The  percent¬ 
age  of  mixing  is  quantified  using  the  following  expression 
[27]: 


%  mixed  =11  — 


x  100 


V(i/AQE^i(Q  -  cf™d)2 

oJmfUiC?  -  c,mixed)2 


(12) 


where  N  is  the  total  number  of  cells,  Q  the  concentration  of 
cell  /,  Cfllxed  the  concentration  of  the  perfectly  mixed  solution 
at  cell  i,  and  Cd  the  concentration  at  cell  i  if  no  mixing  or 

1  _ 1 

diffusion  were  to  occur.  With  an  inlet  velocity  of  0.1  ms  , 
at  the  outlet  of  the  fuel  cell,  Geometries  1-3  resulted  in  14%, 
8%  and  19%  mixing  respectively  (measured  at  the  outlet). 
These  results  indicate  that  particularly  for  Geometry  2,  there 
is  an  opportunity  to  utilize  much  more  fuel. 

The  effect  of  the  inlet  velocity  magnitude  for  Geometry 
2  is  shown  in  Fig.  8.  With  respect  to  an  analysis  of  fuel  uti¬ 
lization  and  mixing,  reducing  the  inlet  velocities  is  effec¬ 
tively  equivalent  to  increasing  the  length  of  the  microfluidic 
fuel  cell.  Decreasing  the  inlet  velocity  from  0. 1  to  0.02  m  s-1 
causes  the  fuel  utilization  to  increase  non-linearly  from  8% 
to  23%.  This  increase  in  fuel  utilization  is  highly  signifi¬ 
cant  in  the  context  of  values  provided  in  previous  studies 
[7,8].  The  mixing  region  at  the  outlet  reached  20  \xm  wide, 
which  is  equivalent  to  the  separation  distance  between  the 
electrodes.  Decreasing  the  inlet  velocity  further  would  result 
in  fuel  cross-over.  It  is  noteworthy,  however,  that  the  per¬ 
centage  mixed  at  the  outlet  increased  only  3%  (from  8%  for 
the  0.1  ms-1  inlet  velocity  to  1 1  %  for  the  0.02  m  s  1  inlet 
velocity),  indicating  that  further  improvements  in  fuel  uti¬ 
lization  are  possible  by  extending  a  tapered  electrode  beyond 


Fig.  7.  Depletion  and  mixing  of  HCOOH  and  O2  in  the  microfluidic  fuel  cell  with  Geometry  3:  (a)  the  location  of  the  X-Z  cross  section  of  interest,  (b)  HCOOH 
mass  fraction  contour  plot  showing  diffusive  mixing  between  the  fuel  and  oxidant  streams  without  electro-oxidation  reactions,  (c)  HCOOH  mass  fraction 
contour  plot  showing  combined  diffusive  mixing  and  electro-oxidation  at  the  anodes  (top  and  bottom,  left),  and  (d)  O2  mass  fraction  contour  plot  showing 
combined  diffusive  mixing  and  reduction  at  the  cathodes  (top  and  bottom,  right). 
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Fig.  8.  Effect  of  increasing  the  inlet  velocity  on  fuel  utilization  and  diffusive  mixing  region  width  for  Geometry  2,  measured  at  the  outlet  of  a  6  mm  long 
microfluidic  fuel  cell. 
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Fig.  9.  A  Schematic  diagram  of  the  microfluidic  fuel  cell  system  with  Geometry  2  extended  with  tapered  electrodes. 


the  6  mm  length  (while  maintaining  channel  cross-sectional 
dimensions). 

The  schematic  of  Geometry  2  with  extended  tapered  elec¬ 
trodes  is  given  in  Fig.  9.  The  concept  is  to  taper  the  electrodes 
to  match  the  growth  of  the  mixing  region,  and  thereby  exploit 
unused  and  unmixed  fuel  while  avoiding  fuel  cross-over.  With 
an  inlet  velocity  of  0.02  ms-1,  this  electrode  design  results 
in  a  fuel  utilization  increase  from  23%  to  52%.  Fig.  10  il¬ 


lustrates  the  formic  acid  and  oxygen  concentration  boundary 
layers  developed  at  the  electrode  surfaces,  along  with  the 
diffusive  mixing  regions  for  the  extended  electrode  geom¬ 
etry.  Compared  to  the  species  mass  fraction  contour  plots 
shown  in  Fig.  6,  where  the  fuel  cell  length  extended  to  only 
6  mm,  the  contour  plots  shown  in  Fig.  10  illustrate  that  the 
fuel  utilization  has  been  significantly  improved  with  the  ex¬ 
tended  geometry  and  tapered  electrodes.  Fig.  1 1  illustrates 


Fig.  10.  Depletion  and  mixing  of  HCOOH  and  O2  in  the  microfluidic  fuel  cell  with  Geometry  2  extended  with  tapered  electrodes:  (a)  the  location  of  the  X-Z 
cross  section  of  interest,  (b)  HCOOH  mass  fraction  contour  plot  showing  diffusive  mixing  between  the  fuel  and  oxidant  streams  without  electro-oxidation 
reactions,  (c)  HCOOH  mass  fraction  contour  plot  showing  combined  diffusive  mixing  and  electro-oxidation  at  the  anodes  (top  and  bottom,  left),  and  (d)  O2 
mass  fraction  contour  plot  showing  combined  diffusive  mixing  and  reduction  at  the  cathodes  (top  and  bottom,  right). 
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Fig.  11.  Depletion  of  HCOOH  in  the  microfluidic  fuel  cell  with  Geometry  2  with  tailored  electrodes:  (a)  tailored  electrode  placement,  (b)  HCOOH  mass 
fraction  contour  plot  on  the  X-Y plane  at  z  =  13.75  |xm  (quarter  plane),  (c)  HCOOH  mass  fraction  contour  plot  on  the  X-Y plane  at  z  =  27.5  |xm  (mid-plane)  and 
at  X-Z  planes  along  the  downstream  channel  (y  =  2,  6,  12  and  18  mm).  Length  scales  were  adjusted  to  illustrate  the  depletion  of  HCOOH  in  the  downstream 
direction. 


the  depletion  of  formic  acid  as  the  fuel  flows  downstream. 
Fig.  1 1(a)  shows  the  tailored  electrode  placement.  Fig.  1 1(b) 
is  the  formic  acid  mass  fraction  contour  plot  on  the  X-Y  plane 
at  a  depth  of  z  —  13.75  pan  (quarter  plane).  Fig.  11(c)  is  the 
formic  acid  mass  fraction  contour  plot  on  the  X-Y  plane  at  a 
depth  of  z  =  27.5  jam  (mid-plane).  Due  to  the  concentration 
boundary  layers  near  the  top  and  bottom  walls  of  the  channel, 
the  formic  acid  concentration  is  more  depleted  in  the  quar¬ 
ter  plane  compared  to  the  mid-plane.  Fig.  11(c)  also  shows 
formic  acid  mass  fraction  contour  plots  at  X-Z  planes  along 
the  downstream  channel  at  positions  y  =  2,  6,  12,  and  18  mm. 
Similar  to  the  fuel  mass  fraction  contour  in  the  X-Y  plane, 
the  contours  in  the  X-Z  planes  also  illustrate  the  depletion  of 
formic  acid  as  the  fuel  travels  downstream. 

Geometry  2  yields  a  total  current  of  0.059  mA,  and  when 
this  geometry  is  extended  with  tapered  electrodes  the  total 
current  increases  to  0.13  mA  corresponding  to  an  average 
current  density  of  1.1  mA  cm-2.  It  is  noteworthy  that  the  ex¬ 
perimentally  measured  output  power  would  be  expected  to 
be  less  than  the  theoretical  prediction  due  to  the  effects  asso¬ 
ciated  with  non-constant  overpotentials  and  fluid  properties. 
Furthermore,  the  resistivity  of  the  fuel  cell  would  change  as 
the  liquid  flows  downstream  where  the  electrolyte  concentra¬ 
tion  weakens.  These  results  indicate  that  extending  the  length 
of  the  fuel  cell  geometry  with  tapered  electrodes  provides  a 
larger  total  current  and  an  improved  use  of  reactants. 


5.  Conclusions 

This  work  is  one  of  the  first  CFD  analyses  for  this  emerg¬ 
ing  technology.  Recent  experimental  studies  have  estab¬ 
lished  proof-of-concept  results;  however,  poor  fuel  utilization 
greatly  limits  the  performance  of  current  designs.  Improving 
fuel  utilization  while  minimizing  fuel-oxidant  mixing  in  the 
microfluidic  fuel  cell  was  the  focus  of  this  study.  A  concise 
electrochemical  model  of  the  key  reactions  and  appropriate 
boundary  conditions  for  the  microfluidic  fuel  cell  have  been 
developed  and  implemented  into  the  numerical  model.  A  nu¬ 
merical  investigation  of  the  coupled  flow,  species  transport 
and  chemical  aspects  provided  insight  into  the  running  pa¬ 
rameters  and  cell  geometry  required  to  achieve  significantly 
improved  performance. 

The  geometry  of  the  cell  plays  a  critically  important  role 
in  the  performance  of  the  cell,  which  was  limited  by  the  mass 
transport  of  reactants  through  the  concentration  boundary 
layers  to  the  electrodes.  Three  geometries  were  investigated. 
The  straight  and  rectangular  channel  geometries  simulated 
in  this  study  are  in  keeping  with  current  designs.  The  rectan¬ 
gular  geometry  with  a  high  aspect  ratio  in  the  cross- stream 
direction  is  the  most  promising  design  for  the  microfluidic 
fuel  cell,  as  the  reactants  demonstrate  the  least  percentage  of 
mixing  and  the  best  fuel  utilization  at  the  outlet  of  the  fuel 
cell.  With  the  same  hydraulic  diameter  as  the  square  geome- 
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try,  the  overall  fuel  utilization  increases  by  more  than  a  factor 
of  two  over  the  square  geometry. 

Lowering  the  inlet  velocity  results  in  a  significant  im¬ 
provement  in  fuel  utilization  with  a  minimal  increase  in  cross¬ 
stream  mixing  at  the  outlet.  Decreasing  the  inlet  velocity  from 
0.1  to  0.02 ms-1  causes  the  fuel  utilization  to  increase  non- 
linearly  from  8%  to  23%,  and  only  causes  an  increase  of  3% 
in  cross-stream  mixing  at  the  outlet.  Power  output  require¬ 
ments  and  the  diffusive  mixing  of  fuel  and  oxidant,  how¬ 
ever,  ultimately  limit  the  degree  to  which  the  velocity  can  be 
reduced.  Based  on  our  analysis,  a  tapered-electrode  design 
was  suggested.  Numerical  simulations  indicate  that  signifi¬ 
cant  improvements  result  from  the  tapering  of  the  electrodes 
towards  the  outlet  of  the  microfluidic  fuel  cell,  while  main¬ 
taining  constant  channel  cross-sectional  dimensions.  By  tai¬ 
loring  the  electrode  shape  to  the  reactant  concentration  profile 
for  an  extended  cell  geometry,  the  fuel  utilization  increased 
from  23%  to  52%,  producing  an  average  current  density  of 
1.1  mA  cm-2. 
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